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Abstract 

Both dynamic and static approaches are proposed and investigated for controlhng the optical 
phase of a p-polarized Ught wave that is guided through a surface-patterned metalhc structure with 
subwavclength features. For dynamic control, the important role of photo-excited electrons in a 
slit-embedded atomic system with field-induced transparency (FIT) is discovered within a narrow 
frequency window for modulating the intensity of focused transmitted light in the near-field region. 
Based on the electromagnetic coupling. This is facilitated by surface plasmons between the two 
FIT-atom embedded slits. The near-field distribution may be adjusted by employing a symmetric 
(or asymmetric) slit configuration as well as by either a small or large slit separation. Furthermore, 
the cross-transmission of a light beam is predicted as a result of strong coupling between optical 
transitions within embedded FIT atoms and surface plasmons. For static control, the role of 
surface curvature is obtained for focused transmitted light passing through a Gaussian-shaped 
metallic microlens embedded with a linear array of slits, in addition to a negative light-refraction 
pattern, which is associated with higher-diffraction modes of light, under a large angle of incidence 
in the near-field region. Most interesting, however, this anomalous negative light-refraction pattern 
becomes suppressible with leaked higher-order waveguide modes of light passing through a very 
thin film. At the same time, it is also suppressible with a reinforced reflection at the left foothill 
of a Gaussian-shaped slit array for the forward-propagating surface-plasmon wave at large angles 
of incidence. A prediction is given for near-field focusing of light with its sharpness dynamically 
controlled by the frequency of the light in a very narrow window. Moreover, a different scheme 
based on Green's second integral identity is proposed for overcoming a difficulty in calculating the 
near-field distribution very close to a surface by means of a finite-difference-time-domain method. 
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I. INTRODUCTION 



Dipole-coupling based optical transition of electrons in an atom is usually very sharp and 
strong under a resonant frequency of the light. As a result, the coherent electron dynamics 
occurs only within an extremely-narrow frequency window. If a laser is used to resonantly 
excite a multi-level atomic system, the quantum interference will simultaneously occur in the 
optical absorption spectrum as well as in the spontaneous emission spectrum at the same 
time. A quantum electrodynamic theory, that treats both photons and electrons in second 
quantization, is compulsory for looking into the radiative decay of excited electrons in an 
atomic system^l^. In addition to the well-known diagonal radiative-decay time of atoms, there 
exists another less-known off-diagonal radiative- decay coupling (ODRDC) that becomes sig- 
nificant as two or more electronic transition energies are very close.lZElExphcitly, the ODRDC 
phenomenon is associated with a nearly-resonant reabsorption of a spontaneously emitted 
photon from one radiative decay process by another upward stimulated optical transition of 
electrons. Consequently, any two nearly-degenerate electronic transitions in the presence of 
ODRDC will be coupled and interfere with each other, leading to a zero absorption linel^ 
in the spectral region where two absorption peaks are overlapped. 

The collective excitation of electrons in metals can couple to light to form surface-plasmon 
polaritons through dielectric-function dependent boundary conditions for electromagnetic 
fields on both sides of an interface. The incoherent dynamics of excited surface-plasmon 
polaritons usually shows up in a broad frequency range owing to the high optical loss of 
metals. In recent years, surface-plasmon polaritons'^ and local surface plasmons'^ as well, 
which are localized at a dielectric-metal interface, have acquired a great deal of attention'^. 
The experimentally observed extraordinarily-high transmission of a p-polarized light beam 
passing through a two-dimensional hole array embedded on a metal film in the subwave- 
length regimell^^ has been proved to depend on both the array period and film thickness 
in a sensitive way'^^'^. In addition, studies on the surface-plasmon mediated light trans- 
mission by a subwavelength structure on a metal surface have also attracted a great deal of 
attention!™!. If 

a single slit is perforated on an optically-opaque metal film deposited on a 
dielectric substrate, it is discovered that the induced surface plasmons on the front side of 
the film will be coupled to the backside by intra-slit interferences'^'^. It is also found that 
not only the intra-slit interference but also the inter-slit interference in a slit array is a vital 
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component for the transmission of the excited surface plasmons passing through the shts'^. 
Furthermore, the inter-sht interference is shown to relate only to the surface-wave nature 
and not necessarily to the surface-plasmon wave, because this interference survives even for 
a perfect electric conductor. 

It is known that surface patterning of a metal film at the nanoscale, e.g., slits, holes, and 
even surface corrugations, can provide flexible control of light incident on it, for example, 
fiat metallic lenses have been constructed with nanoslits perforated on a thin metal film to 
demonstrate beam deflection and focusing based on specific phase retardation controlled by 
modulated slit depths'^, slit widths slit separations-^^ and slit dielectric filling'I^. More- 
over, surface plasmons have been identified as the essential ingredient in these phenomena, 
and the Fabry-Perot resonance of surface plasmons contributes greatly to the enhanced 
transmission and localized electric-field in the slit perforated metallic film'^^ES. The config- 
urations of a single aperture smaller than the incident wavelength and an array of nanoslits 
perforated on a metallic film have attracted a great amount of attention due not only to 
extraordinary optical transmission of light^, but also to their beaming effect in the far 
gg}(j[25H29] g^j^^ ^Yie "lens-like" properties'^^'^ for focusing normally-incident light in the inter- 
mediate region within a certain frequency range. In the near-field region, the propagation 
of a light wave guided through a metallic structure in the subwavelength regime is deter- 
mined by light diffraction, instead of by Snell's law for a light ray. The optical phase of 
a guided light wave is a product of the host-material refractive index and the geometrical 
distance traveled by the light inside the host. The frequency-dependent refractive-index 
function reflects the dynamics of photo-excited electrons in the time domain, which can 
be dynamically controlled through an optical resonance between a pair of discrete energy 
levels in addition to a weak Kerr optical nonlinearity, such as embedded atoms or quantum 
dots, within a very narrow frequency window. For this dynamic-control approach, strong 
nonlinear coupling of electrons to an incident light needs to be taken into account, e.g., the 
coupled Maxwell-Bloch equations for including nonlinear optical response of photo-excited 
electrons inside atoms to a spatially-inhomogeneous total electromagnetic field have to be 
solved. On the other hand, the geometrical distance traveled by the light can be statically 
controlled by changing the surface profile of a host material at its interface with either air 
or a substrate, e.g., a Gaussian-shaped metallic microlens embedded with a finite nanoslit 
array. For this static-control approach, only the Maxwell equations with the inclusion of a 
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linear optical response of photo-excited electrons in conductive materials are required. 

The motivation for the current research is stated as follows. First, it is well know that the 
slow and incoherent dynamics of surface-plasmon excitation can produce a very strong near 
field close to a surface and a broad transmission peak in the spectrum. On the other hand, the 
fast and coherent dynamics originated from stimulated transitions of electrons inside atoms 
can lead to a strong absorption peak in which a very narrow transparency window is buried. 
Therefore, the combination of these two quite different dynamics is expected to introduce a 
new mechanism for near-field focusing of light with its sharpness dynamically controlled by 
the frequency of the light in an extremely narrow window in the spectrum. Second, for a 
plane-wave guided through a metallic nanoslit array in the subwavelength regime, it is also 
well known that only the near field can survive as the incident angle of light becomes very 
large. Although a flat lens composed of nanoslits perforated on a thin metal film can be 
applied to focusing normally-incident light, a long tail along the lens axis in the light- focusing 
pattern is unavoidable. The use of a curved nanoslit array embedded on a metal film, on 
the other hand, can provide a very sharp and configuration-dependent focus spot without a 
longitudinal tail even at large angles of incidence. In this sense, the near-field focusing of a 
plane- wave light can be controlled statically mainly through a change of the surface curvature 
of a metallic lens. As far as the near-field distribution very close to a surface is concerned, 
the commonly-employed simple finite-difference-time-domain method drops memory effect 
in the time domain and tends to leave out details in the field-distribution pattern very close 
to the surface, which requires inclusion of very large transverse wave numbers (or very fine 
spatial meshes in the direction parallel to the surface) in numerical calculations. However, 
the analytic calculations based on Green's second integral identity in this work do not suffer 
from this drawback. 

In this paper, we investigate both dynamic and static control of the optical phase for a 
p-polarized light wave guided through a subwavelength surface-patterned metallic structure. 
For dynamic control, we explore the coherent dynamics of photo-excited electrons in a slit- 
embedded atomic system with field-induced transparency (FIT) within a narrow frequency 
window, which is further superposed on the incoherent dynamics of surface plasmons in a 
broad frequency range. For the static control, we investigate the effect of surface curvature on 
focused transmitted light passing through a Gaussian-shaped metallic microlens embedded 
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with a finite array of nanoslits under a very large angle of incidence in the near-field region. 



The rest of the paper is organized as follows. In Sec.[TT| we introduce two model systems, 
respectively, for dynamic and static control of the optical phase of a guided p-polarized light 



wave. In Sec. Ill, numerical examples are presented and discussed in detail for dynamic and 
static optical-phase control in these two model systems. Finally, a brief summary of the 



main results of this paper is given in Sec. IV 



II. MODEL SYSTEMS FOR OPTICAL-PHASE CONTROL 

In this section, we will introduce two model systems associated with the dynamic and 
static control of the optical phases, respectively. Based on these two introduced model 
systems, numerical calculations are performed to demonstrate both the coherent electron 
dynamics in embedded atoms randomly-distributed within a nanoslit and the incoherent 
dynamics associated with interference effects of resonantly excited surface-plasmon waves. 
Also shown is the suppression of interference by reinforced reflection at the "foothill" of a 
curved slit array under large angles of incidence. 



A. Dynamic approach 

By taking into account of the dynamics of photo-excited electrons within embedded 
atoms, the model system introduced for dynamically controlling the optical phase is il- 
lustrated in Fig.jlj The complete theory beyond a so-called "diagonal" approximation'^Dl'^, 
which is based on the modal expansion and the surface-impedance-boundary condition, 
for calculating the spatial distribution of the total electromagnetic field as well as the 
transmission and reflection spectra in such a model system has already been reported 
in Ref. [18]. Here, the dielectric constant kj for the nanoslit-filled material in Ref. [18] 
should be replaced by a frequency-dependent effective one with the inclusion of randomly- 
distributed atoms (random in position but with the same locally-spherical symmetry) em- 
bedded within the nanoslit-filled dielectric medium. Explicitly, we require the substitution 
Kj — )■ ej{ijj) = Kj + a-[{uj) in Ref. [H] with the Lorentz ratio of photo-excited electrons given 
byia 
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aiiuj) 



(1) 



where Cj is the concentration of embedded atoms in the jth nanosht, hu is the incident pho- 
ton energy, = 2d^^EQ/h for z/ 7^ /i = 1, 2, 3 (three-level system) is the Rabi frequency, 
dj,^ is the dipole moment for the optical transition of electrons from a lower energy level 
|z/ > to a higher energy level |/i >, 
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, and = ^ / d\^Eir\^) . (2) 

Oj JSj 

In Eq. (j2]), the spatial averages of the product £'(r||) 0(r||) and E{y\\) itself are performed 
within the jth-nanoslit with an area iSj, y\\ is a two-dimensional position vector within the 
xz-plane, E{y\\) is the magnitude of the electric-field component, and 0(r||) represents its 
structure factor in the xz-plane. 

In order to relate the electric-field vector E(r||) to the solution of the Maxwell equations 
in Ref. [TH], we consider a p polarization for the incident light as an example for this paper. 
Therefore, the total magnetic-field vector H(r||) in the direction perpendicular to the xz- 
plane, as shown in Fig.[l| can be simply written as H(r||) = [0, Hy{x,z), 0]. This gives rise 
to an expression for the electric-field vector 



Efrii) = 



as well as its structure factor, given by. 



dHy{x,z) ^ dHy{x, z) 
(9z ' ' dx 



0(rii) = 



(3) 



1 dHy{x,z)/dx — dHy{x,z)/dz 
V3 ^J[dHy{x, z)/dxf + [dHy{x, z)/dzf 



Finally, the density-matrix elements pj^^iu)) for the jth-nanoslit, which are introduced in 
Eqg for /i, z/ = 1, 2, 3, can be calculated by assuming a three-level atomic system with 
the property of field-induced transparency^^, as shown in Fig.[2]^a). The details, which are 
based on the optical Bloch equations with quantum interference from the ODRDC, for cal- 
culating the density-matrix elements pj^i,{uj) in such a system have already been formulated 
in Refs. |[3[|6j. The iterations are required for solving the coupled nonlinear Maxwell-Bloch 
equations and getting the self-consistent solutions with respect to Hy{y,z) and p^j,(w) at 
the same time. 
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B. Static approach 



By including a nontrivial surface profile to the interface of a host material, the model 
system introduced for statically controlling the optical phase is illustrated in Fig. [6] The 
in-detail formalism, which is based on Green's second integral identity in the x2;-plane'^, for 
calculating the spatial distribution of the total electromagnetic field as well as the transmis- 
sion and reflection spectra in such a model system has already been reported in Refs. [32|33] . 
The top surface profile in Ref. j33] should be replaced by 

[ A exp (-^) E exp [-^^1 , for \x\ < D,/2 

[ Qo , for |a;| > Do/2 

where j is the nanoslit index, d is the array period, 2b represents the width of nanoslit- 
separation layer, 2N + 1 is the total number of nanoslits with width ~ (rf — 2b). In addition, 
Dq is the microlens aperture size, Qq + a ^ Qo is the metal film thickness, (a is much 
smaller than the incident wavelength Aq for strong photon tunneling), Ao and 2ao are the 
amplitude and width for the nanoslit-depth modulation. Additionally, the bottom surface 
profile in Ref. |33] is replaced by a flat one, i.e., ^2{x) = —a. For |a;| <^ 1, the quadratic 
term in the expansion of ^i(x) in Eq. Q leads to an effective curvature-based focal distance 
fo ~ 0-0/(2 A) measured from the exit flat surface. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

In our numerical calculations, the metal is chosen to be silver (gold) for dynamic (static) 
control, and its dielectric functions eui^jj) with different frequency u of the light are taken 
from the work by Johnson and Christy'^. 



A. Quantum-interference effect in dynamic control 

For a filled metallic nanoslit array shown in Fig.[T} the parameters taken for the numerical 
calculations in Figs. [2]-|5] are: eL = 1 for air, eR = 2.13 for glass, 2d = 1.155 /xm for film 
thickness, and = 0° for normal incidence. For the stand-alone three-level atomic system, 
we take: /32i,i2M2 = /33i,i2M2 = /32i,i3M2 = /33i,i2M2 = 0.1, 70M2 = 0.1, d23 = 0, 
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1^12 = c?i3, S = {ujp~ uj2i)/uj32, ^vi/^?.2 = ^^12^2 = ^p/(^32, = 2di2Ep/h with Ep being 
the probe-field amplitude. Furthermore, for the embedded atoms within the nanoslit-fiUed 
dielectric material, we choose the same atom concentration Cj = 1 x 10^^ cm~^ for all the 
nanoslits, hujs2 = 0.1 meV, and di2/e = dis/e = 1 A. Other parameters as well as the changes 
of these parameters will be given in the figure captions. For all plots in Figs.|3]-|5j the color 
scale is kept the same from (blue) to 1 (red). 

Figure [2] presents the scaled Lorentz-ratio function af^{u}p) (2eoc<;32-E'p/cat^^^p) for a stand- 
alone atomic system with field-induced transparency, as shown in Fig.|2](a), as a function 
of the dimensionless detuning parameter 6 = {up — a;2i)/w32 for different strengths of a 
probe field Ep, given by ^lp/ujs2 = 0.067 [weak, in Fig.[2](b)], 0.2 [medium, in Fig.|2](c)] and 
0.67 [strong, in Fig.|2](d)]. When Ep is weak, we find a zero absorption from Im[a;L*(ci;p)] 
in Fig.[2]^b) at S = 0.5, i.e., the field-induced transparency (FIT) owing to quantum inter- 
ference between direct and an indirect optical-transition phase based on ODRDC'^. This 
transparency is accompanied by two absorption peaks at both sides oi 6 = 0.5. In addition, 
corresponding changes in Re[a;L*(co'p)] with respect to these two absorption peaks are seen 
with a very large negative slope for d{Iie[a^ (up)]} / dup within a very narrow transparency 
window. As Ep is gradually increased from Fig.|2]^c) to Fig.|2]^d), the strength of these two 
absorption peaks greatly decreases while the width of them (power broadening) significantly 
increases at the same time. The peak reduction will continue with increasing Ep until the 
complete suppression of absorption peaks is achieved at a saturation limit. The results in this 
figure clearly demonstrate a rapid dynamic change of a^{up) and a very strong probe-field 
nonlinear ity within an extremely-narrow frequency window Aup ~ 0032 ^ 7o- 

When the FIT atomic system in Fig.[2]^a) is embedded into a single nanoslit, the probe 
field will be replaced by a guided light wave. At the same time, the electron dynamics 
within a narrow-frequency window for an isolated atomic system is also inserted into a 
relatively broad frequency range for the dynamics of surface-plasmon polaritons. Moreover, 
the FIT atoms will be coupled to surface-plasmon modes in a nanoslit. For a very narrow 
and long nanoslit in the deep-subwavelength regime, we know for p-polarization incidence 
only the lowest waveguide mode can be transmitted through the nanoslit. We display in 
Fig. [3] the self-consistently calculated transmission coefficient in the upper-left panel and 
color maps for \Hy{x,z)\'^ in the other three panels with various detuning parameters 6 (or 
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equivalently, different FIT atoms) in tlie rest tliree panels below the saturation limit. A high 
light reflection by the front surface of a silver film, as well as the induced surface-plasmon 
propagation along this surface, can be seen from Fig.|3]^b) — (d), but the pattern of \Hy{x, z)p 
within the nanoslit remains the same except for intensity variation. It is important to 
note that the nanoslit has been prepared for one of the passing states satisfying 8d = 
|2m— 1 1 Ao / y/Ko with m = 0, ±1, ±2, ■ ■ ■ resulted from dual-wave constructive interference'^. 
Compared with the medium {6 = 0.25 or ijJ2i = 27tc/Xq — c<;32/4) and with the strongest 
{6 = 0.5 or u)2i = 27rc/Ao — CJ32/2) light transmissions in the lower-left and lower-right panels 
of the figure, the transmitted light at 5 = (or U21 = 27tc/Xq) in the upper-right panel of 
the figure is the weakest owing to maximized absorption at this particular 6 value, as can be 
verified by the calculated transmission coefficient in the upper-left panel of this figure. In this 
sense, the intensity of transmitted light can be self-modulated (or dynamically controlled) 
in a narrow-frequency window by embedding atoms within a nanoslit of a silver film, and 
the modulation effect can be enhanced by increasing the concentration of embedded atoms 
in the nanoslit. 

Figure [4] displays the calculated color maps for \Hy{x,z)\'^ with double nanoslits in four 
different cases. In the upper-left and upper-right panels of this figure, we consider a symmet- 
ric configuration of double nanoslits with 6 = and 6 = 0.5, respectively, for both nanoslits. 
Similar to the color maps in Fig. [3] for a single nanoslit, we observe a small reduction in 
transmitted light through nanoslits for 5 = due to a relatively large absorption. Because 
of the electromagnetic coupling between two nanoslits, a hot spot occurs in the reflection 
side at the front surface between two close-by nanoslits. Moreover, two candle-fiame-like 
spots at the nanoslit exits are significantly deformed in shape and bent down from the mid- 
dle point of an individual nanoslit towards the middle point between two nanoslits, which 
is quite different from that seen in Fig.[3j In the lower- left panel of the figure, we consider 
an asymmetric configuration for double nanoslits with 6 = and 6 = 0.5 for the lower and 
upper nanoslits, respectively. Although the hot spot still exists at the front surface for the 
reflection because of established electromagnetic coupling between two different nanoslits, 
the bent-down candle-flame-like spot at the upper-nanoslit exit becomes brighter than that 
at the lower-nanoslit exit for this asymmetric case. Furthermore, when these two asym- 
metric nanoslits are moved apart from each other, as shown in the lower-right panel of the 
figure, the electromagnetic coupling between nanoslits no longer exists, which leads to a dis- 

10 



appearance of the hot spot in the reflection and a recovery of two deformed candle-flame-hke 
spots. The lack of couphng between nanoshts also leads to the enhanced surface-plasmon 
propagation along the back surface of the silver film resulted from the resonance between 
the incident light and the surface plasmon at the metal-glass interfacel^. 

The results in Figs. [3] and |4] correspond to a simultaneous illumination of all nanoslits 
by a p-polarized light in the plane- wave form. On the other hand, we find from Fig. [5] 
that the nonlocal dynamic effects from embedded atoms in an unilluminated nanoslit, as 
a feedback, can affect the transmission through a separated illuminated nanoslit without 
embedded atoms. Here, the bottom one of the two lower non-embedded nanoslits is illumi- 
nated by a narrow Gaussian light beam while the top unilluminated nanoslit is embedded 
with FIT atoms. At Aq = 620 nm for the two upper panels of this figure, the bottom il- 
luminated nanoslit stays in an intermediate state'^ having a medium direct transmission 
of light and a stiff local reflection, while the top unilluminated and fllled nanoslit still ac- 
quires a relatively high cross-transmission coefficient much larger than that of the middle 
unilluminated nanoslit. Moreover, the stiff reflection enhances the propagation of a resonant 
surface-plasmon mode at the metal-air (front) interface at this wavelength. From the upper 
left panel of the flgure, it is important to mention that the electron dynamics in the top 
nanoslit has been mapped onto the dynamics for the direct transmission of light through 
the bottom non-embedded nanoslit, where two Fabry-Perot-type peaks right above two dips 
for the top nanoslit are observed with their replicas in the bottom nanoslit due to inter-slit 
coupling. This coupling is facilitated by the resonant surface-plasmon mode localized at 
the metal-air interface. This self-induced transmission for the embedded nanoslit is found 
to be related to the sole- wave interference'^ between the entry and exit edges of a nanoslit 
with a finite value for the average refiection coefficient (a large finesse increasing with the 
film thickness 2d), where the requirement for the Fabry-Perot constructive interference is 
fulfilled by a sharp increase in Re[aL(i^)] right above 5 = or 5 = 1. After Aq is switched 
to a longer wavelength at 990 nm, as shown in the two lower panels of this figure, the bot- 
tom illuminated nanoslit changes from an intermediate state to a passing state with a large 
direct transmission of light and a softened local refiection. In this case, the propagation of 
the previous resonant surface-plasmon mode at the metal-air (front) interface is suppressed 
at this longer wavelength, while a new resonant surface-plasmon mode at the metal-glass 
(back) interface is induced. As a result, two Fabry-Perot-type peaks for the top nanoslit 
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become invisible owing to the cutoff of surface-plasmon propagation at the front interface 
and their rephcas in the bottom nanosht are also greatly weakened. 



B. Curvature effect in static control 



For a shaped metallic microlens embedded with a finite nanoslit array, as shown in Fig. [6], 
the parameters used for the numerical calculations in Figs. [7]- 10 are taken to be: = 14, 



d = 60 nm, 26 = 6 nm, = 400 nm, 2aQ = Dq = 1.8 fim, Qq = 152 nm, a = 4 nm, egub = 13, 
and eair = 1- Various parameters are given in the figure figure captions. 

By neglecting the small surface-plasmon-polariton phase, the spot-focusing condition to 
the leading-order requires that the difference in the optical phases corresponding to any 
pair of nanoslit channels should be equal to a multiple of the incident wavelength Aq'^^^. 
Therefore, for a given arbitrary shape function g{x) assigned to the front surface, it is found 
to satisfy the following constraint under normal incidence of light 

«M= , ,f^' ^ (v^-/.)- , ^ . (5) 



-air 



where x = jd for j = 0, ±1, ±2, ■ • ■ , ±A^, m = 0, ±1, ±2, ■ ■ ■ is an integer, /q repre- 
sents the designed short distance of a focused spot measured from the flat back surface. 



yeefr(Ao) = /3siit(Ao)/fco stands for the effective complex dielectric function of a nanoslit 
waveguide including the surface-plasmon contribution, ko = 27r/Ao, and the complex 
wavenumber /3siit(Ao) can be determined from the nanoslit-waveguide boundary condition 
for a p-polarized plane wave. By choosing g{x) = ~ ^li^) ~ -^o^^/'^o \^\ ^ '^O; "we 
find from Eq. ([s]) under m = 0, to the leading order, a nanoslit-independent focal distance 



X J~n(X} X 

JO '- 



2To{x) 2J^o{x) 



ecfr(Ao 



(6) 



2/C, 



where J'oix) = [(\/eeff(Ao) - ^/e^T) / ^/^]{K^oX^ / ctq) , = A/o-q, and a larger Kq value 
corresponds to a small value for the focal distance /q. 

Served as a purpose for comparison and a starting point for the static-control approach, we 
display in Fig.Hthe calculated color maps of \Hy{x, z)]"^ for an unshaped gold film embedded 



with a finite nanoslit array at four different angles of incidence. In order to obtain a sharply- 
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focusing spot for the transmitted light at a very short distance /o in such a case, the nanosht- 
center positions Xj with j = ±1, ±2, ■ ■ ■ (and Xq = 0) for a p-polarized normal incidence 
must satisfy a relation, i.e., x'j = (|j|Ao/y'ecfr(Ao) )^ + 2/o|j|Ao/y'eeff(Ao). For a periodic 
nanoslit array with xj = jd and j = 0, ±1, ±2, • • • , ±A^, on the other hand, we expect an 
elongated focusing spot in the z direction, as seen in the upper-left panel of the figure at 
6i = 0". Additionally, the elongated spot outside of the near-field region is associated with 
the specular (m = 0) mode, while those symmetric side robes with respect to Xq = are 
resulted from m = ±1, ±2, • • • higher-order diffraction modes within the near-field region. 
When 6i = 30° in the upper-right panel of the figure, the specular mode is tilted accordingly 
similar to Snell's law. At the same time, the spatial distribution of higher-order diffraction 
modes change and become asymmetric with suppressed diffraction modes at the right side 
for m = 1, 2, ■ ■ - . As 6'j is further increased to 45° in the lower-left panel of the figure, the 
specular mode is tilted even more, as expected, and a series of hot spots in the near-field 
region is pulled back to the exit surface and squeezed towards the right edge of the nanoslit 
array at the same time. Eventually, the specular mode becomes completely suppressed 
outside the near-field region in the lower-right panel of the figure at 6i = 60°, leading to a 
negative light-refraction pattern supported by three hot spots corresponding to one specular 
mode and two diffraction modes with m = — 1 and —2 in the near-field region. 

For a shaped gold-film microlens embedded with a nanoslit array, we present in Fig. [8] 
the color maps of \Hy{x,z)\'^ with various 6i values. Compared with the results in Fig.jT], 
the elongated focusing spot originated from the specular mode in the upper-left panel of 
this figure is greatly shrunk in size and enhanced in strength simultaneously. This is further 
accompanied by reduced intensities for diffraction modes with m = ±1 and ±2. Here, the 
enhancement in focusing capability is attributed to the added compensation in the optical 
phase provided by the geometrical effect from a finite surface curvature for the shaped 
microlens. The observed focal distance (~ 2.5 fim) at 6i = 0° is qualitatively consistent with 



the prediction (/o = 2.8 yum) from Eq. (6) (by simply taking ^ecfT(Ao) ^ iy/e^+ y/e^)/2 = 
2.3 for the lowest p-polarized waveguide mode). The deviation from an ideal spot shape in 
the upper-left panel of the figure comes from the nanoslit contribution in the outer region of 
\x\ > o"o = 0.9 yum. The focusing quality is expected to be improved if only the central part 
of the nanoslit array is under an illumination by a narrow Gaussian light beam instead of a 
wide plane wave illuminating the entire array. The hot specular mode outside the near-field 
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region is greatly cooled down when 6i > 30°, as can be verified by the comparison of the rest 
of the panels of this figure with their corresponding ones in Fig.[7j Moreover, the negative 
light-refraction pattern in the near-field region at 6i = 60°, as shown in the lower-right panel 
of the figure, becomes much cleaner in vision after a weakly focusing spot from the m = —3 
diffraction mode enters into the game. 

As predicted by Eq. ([6]), the focal distance will depend on the surface curvature of a 
shaped microlens for the normal incidence of a p-polarized light. We display the color maps 
of \Hy{x,z)\'^ in Fig.jsjfor a reduced surface curvature (proportional to /Co) as well as for 
a thinner gold film under different angles of incidence. First, the crossover from a flat 
film to a shaped microlens can be understood by a comparison of the results in this figure 
with those in Figs. [7] and [8} including the evolution of an elongated focusing spot from the 
specular mode outside the near-field region in the z direction, the rapid development of side 
robes from the high-order diffraction modes in the near-field region, as well as the change 
in the focal distance. Secondly, it is interesting to note the varied pattern of the near-field 
distribution from the upper-left-panel of the figure and a new tiny focusing spot as well at 
a much shorter distance from the exit surface due to leaked higher-order waveguide modes 
in this case. Finally, the reduced gold-film thickness makes it possible for the contribution 
to light transmission from the leakage of lossy higher-order waveguide modes, which also 
show up in the coexistence of both the positive and negative light-refraction patterns in the 
near-field region at 6i = 60° in the lower- right panel of the figure. Compared with Fig.[7| we 
conclude that the focusing quality can be further improved if lossy higher-order waveguide 
modes can be fully filtered out by employing narrower nanoslits. 

In addition to the surface curvature. Equation ^ also indicates that the focal distance 
will depend on the incident wavelength Aq through the A/eefr(Ao) factor. Figure 
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shows the 



color maps of \Hy{x, z)\'^ aX a resonant surface-plasmon wavelength (air-metal interface) at 
Aq = 630 nm with four values of 6i. For 6i = 0° in the upper-left panel of this figure, we find 
a reduced focal distance /o, which is in agreement with the prediction /o ~ y^Wi/Ao from 
our previous study for a fiat film embedded with a nanoslit array'^. The strong coupling 
of the resonant surface-plasmon modes at two parts of the metal-air interface outside the 
microlens aperature to the nanoslit waveguide modes produces intensive reflections at two 
edges of the array in opposite directions. Moreover, the pattern of the near-field distribution 



14 



for this case is significantly modified by a stronger, smaller and closer focusing spot resulting 
from the concentrated field at the central part of the array. Here, the change in the near- 
field distribution pattern is closely related to the interference between two resonantly excited 
and counter-propagating surface-plasmon waves under the normal incidence. However, the 
strong array-edge reflections are greatly weakened owing to the decreased interference effect 
from two counter-propagating surface-plasmon waves as 9i increases, as can be seen from the 
other three panels of the figure. At 6i — 45° in the lower-left panel of the figure, a negative 
light-refraction pattern is observed, which is supported by both the specular mode and two 
diffraction modes with m = —2 and —3. Interestingly, after removing the concentrated 
field at the middle part of the array at 9i = 60° in the lower-right panel of the figure, a 
similar field-distribution pattern to that under normal incidence reoccurs because of strong 
reflection of one surface-plasmon wave at the left edge of the array, where the diffraction 
modes with m — —2 and —3 at 9i — 60° play the role of those with m = ±1 at = 0°. 



IV. CONCLUSIONS 



In brief, both dynamic and static approaches have been demonstrated for controlling 
the optical phase of a p-polarized light guided through a metallic film with a patterned 
surface in a deep subwavelength regime. For the dynamic control, the modulation of focused 
transmitted-light intensity has been observed within an extremely narrow frequency window. 
The surface-plasmon based electromagnetic coupling between the two nanoslits embedded 
with FIT atoms has shown additional control over the near-field distribution for symmetric 
and asymmetric configurations as well as for small and large separations between the two 
nanoslits. Moreover, an energy-transferring-process based light cross-transmission has been 
shown for the static control. In addition, the curvature-dependent effect on the position of 
a sharply-focused spot from a transmitted light passing through a Gaussian-shaped metallic 
microlens has been demonstrated at all angles of incidence. 

Physically, a new mechanism for near-field focusing of light has been predicted with 
its sharpness dynamically controlled by the frequency of the light in a very-narrow win- 
dow. Technically, a mitigation scheme based on Green's second integral identity has been 
provided for avoiding drawback in numerically calculating near-field distributions by the 
finite-difference-time-domain method. 
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FIG. 1: (Color online) Schematic illustration for a z direction slit array (blue) which extends in 
the y direction. In our notation, Zj and 2ij denote the center position and the slit width of the j-th 
slit, respectively, with j = 0, ±1, • • • , itA^. The regions on the left- and right-hand sides of the slits 
are indicated as Region I and Region III, respectively, with real dielectric constants eL and epj,. The 
region for the slit array is denoted as Region II, and slits are filled with medium (orange) having 
dielectric constant Kj (real) for j = 0, ±1, • • • , . Randomly-distributed atoms (red dots) are 
embedded inside each slit-filled dielectric medium with a concentration Cj for j = 0, ±1, • • • , itA'^ 
and described by an effective, frequency-dependent and complex dielectric function. The depth 
of slits in the x direction is 2d, and eM('^) represents the dielectric function of the metal film 
containing slits. A Gaussian beam is incident on the slit array from the left-hand side with the 
angle of incidence and at a center position z = zq,. The incident wave number is i/eL /co and /3o 
is the incident wave vector along the z direction. 
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FIG. 2: (Color online) Schematic illustration for the effect due to off-diagonal radiative-decay 
coupling in the bare-atom picture [in (a)] and the scaled aL(^p) (2eof^32-E'p/cat^f^p)5 defined through 
Eq. ([T]), as functions of the detuning 6 = {ujp — 0021) / UJ32 in the other panels (b), (c) and (d), where 
Cat is the atom concentration, huji = ej — £i with 0132 ^ i02i and Ej for j = 1, 2, 3 stand for 
the energies of three labeled states in (a). The notations |1 >, |2 > and |3 > in (a) represent 
different energy states in the bare-atom picture, cOp is the probe-field frequency and di2, dis are 
nonzero optical dipole moments, and (3ij^mn denotes the real part of both the diagonal (/32i,i2 and 
/?3i,i3) and off-diagonal (/32i,i3 and /33i,i2) radiative-decay rates. In (b), (c) and (d) both the real 
(index of refraction) and the imaginary (absorption ) parts of aL(wp) have been presented for the 
pure atomic system in (a) with three different given values for the probe-field amplitude Ep, i.e., 
^p/'^32 = 0.067 (weak), 0.2 (medium) and 0.67 (strong), where ilp = disEp/h, where di3 = di2 
and d23 = are assumed. 
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FIG. 3: (Color online) Transmission coefficient calculated self-consistently as a function of 5 in 
the upper-left panel. Also, shown are color maps for \Hy(x, z)p with a single slit. The slit is ffiled 
with a dielectric material (kq = 2.25) which is embedded with randomly-distributed atoms in the 
rest of the panels with 5 = (upper-right), 0.25 (lower-left) and 0.5 (lower-right), respectively. 
The width and position of the slit are Iq = 0.2 fim and zq = 0. Here, the incident-light wavelength 
is Ao = 990 nm and the p-polarized plane-wave amplitude is fJ-oHo = 0.5 G. 



21 





j 8=(0,0.5) 
I 

i Asymmetric 




FIG. 4: (Color online) Self-consistently determined double-slit color maps for \Hy(x, z)p filled with 
dielectric material = ki = 2.25) which are both embedded with randomly-distributed atoms. 
The upper panels correspond to a symmetry configuration having the same embedded atoms in 
both slits, i.e., 5_i = 6i = (upper left) and 5_i = = 0.5 (upper-right). The lower panels are 
for an asymmetric configuration with different types of embedded atoms in two slits, i.e., 6-i = 
and 5i = 0.5 for the lower and upper slits, respectively. In addition, for the two lower panels, we 
display the results for both coupled (lower-left with a small slit separation) and uncoupled (lower- 
right with a large slit separation) slits. The widths of two slits are l^i = li = 0.2 /im. For the 
upper two and lower-left panels, the positions for these slits are Z-i = —0.25 fim and zi = 0.25 fim, 
while Z-i = — 1 fiui and z\ = \ fxm for the lower-right panel. Here, the incident-light wavelength is 
Ao = 990 nm and the p-polarized plane-wave amplitude is HqHq = 0.5 G. 
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FIG. 5: (Color online) Self-consistent calculations of direct (bottom) and cross (middle and 
top) triple-slit transmissions as functions of 6 in the upper-left (Aq = 620 nm) and lower-left 
(Ao = 990 nm) panels. Also plotted are their corresponding color maps for \Hy{x, z)\'^ with triple 
slits. All slits are filled with the same dielectric material (k_i = kq = ki = 2.25) but only the 
top slit is embedded with randomly-distributed atoms having 5 = 0.5 for two right panels. The 
widths for these three slits are i-i = 0.1 ^m and io = ii = 0.2 /im, while the positions for these 
slits are z_i = — 2.2 ^m, zq = 0, zi = 0.8/xm and zq = zi. In addition, only the lowest slit is 
illuminated with light of a Gaussian beam, the incident-light wavelength is Aq = 620 nm (upper 
panels) or Aq = 990 nm (lower panels), and the peak amplitude of a p-polarized Gaussian-beam is 
HoHq = 0.5 G. 
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FIG. 6: (Color online) Diagram representing a shaped metallic lens. A Gaussian-overlapped profile 
z = ii{x) with an embedded slit array is chosen for the surface on the entry side. A flat surface 
z = i2{x) is assumed for the exit side. A p-polarized plane- wave light with out-of-plane magnetic- 
field component H = (0, Hy, 0) is incident with angle of incidence Oi where the dielectric constant 
is Cair = 1- The substrate dielectric constant is esub = 13 on the exit side. A complex dynamical 
dielectric function eM('^) is employed for metal with optical loss included, where fiuj represents the 
incident photon energy. In addition, the metal-film thickness \s Qq + a ^ Qq {a <^ Aq with Aq 
being the incident wavelength), and Dq represents the lens aperture size. 
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FIG. 7: (Color online) Color plots of \Hy(x, z)\'^ for transmitted light through a flat and slit- 
embedded metallic film when Aq = 540 nm. Various angles 9i of incidence are chosen. The color 
scales are from (blue) to 4.5 (red) for Oi = 0°, 30°, to 5 for Oi = 45° and to 4 for Oi = 60°. 
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FIG. 8: (Color online) Density plots of \Hy{x, z)\'^ for transmitted light through a shaped and 
slit-embedded metallic film when Aq = 540 nm for various angles of incidence 9i. The color scale 
ranges from to 5 for 9i = 0°, 60°, to 6 for Oi = 30° and to 7 for for Oi = 45°. 
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FIG. 9: (Color online) Density plots of \Hy{x, z)\'^ for light transmitted through a shaped and 
slit-embedded metallic film. We chose Aq = 540 nm with the half values for and Qq used in 
Fig. [7] and various angles of incidence 6i. The color scale ranges from to 4 for 9i = 0°, 30°, to 
5 for 9i = 45° and to 3 for for 9^ = 60°. 
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FIG. 10: (Color online) Density plots of \iiy{x^ z)p for transmitted light through a shaped and 
slit-embedded metallic film when Aq = 630 nm for various angles of incidence Qi . The color scale is 
from (blue) to 5 (red) for 6*^ = 0°, 60°, to 6 for Oi = 30° and to 7 for for Oi = 45°. 
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